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Abstract 

Metapopulation theory for a long time has assumed dispersal to be symmetric, 
i.e. patches are connected through migrants dispersing bi-directionally without a 
preferred direction. However, for natural populations symmetry is often broken, 
e.g. for species in the marine environment dispersing through the transport of 
pelagic larvae with ocean currents. The few recent studies of asymmetric disper- 
sal concluded, that asymmetry has a distinct negative impact on the persistence 
of metapopulations. Detailed analysis however revealed, that these previous 
studies might have been unable to properly disentangle the effect of symmetry 
from other potentially confounding properties of dispersal patterns. We resolve 
this issue by systematically investigating the symmetry of dispersal patterns 
and its impact on metapopulation persistence. Our main analysis based on a 
metapopulation model equivalent to previous studies but now applied on regular 
dispersal patterns aims to isolate the effect of dispersal symmetry on metapopu- 
lation persistence. Our results suggest, that asymmetry in itself does not imply 
negative effects on metapopulation persistence. For this reason we recommend 
to investigate it in connection with other properties of dispersal instead of in 
isolation. 

Keywords: Connectivity matrix, dispersal network, symmetry, 
metapopulation viability. 



1. Introduction 

Many species are structured in space with dispersal and migration connectin g 
local populations into metapopulations (|Levinsl llfffifj lHanski and Gilpinl . ll997h . 



The fundamental dynamics of metapopulations are determined by local extinc- 
tion, dispersal from the local populations, and colonisation success leading to the 
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establishment of new sub-populations ( Levins , 1969() . Metapopulation dynamics 
may dete rmine a range of ecological and evolutionary aspects including popula- 



tion size ( Gyllenberg and Hanski . 1992 ). persist ence (Rov et al. . 20051) . spatial 
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distribution (IRov et all 120081 ) , epidemic spread (IMcCallum and Dobsonl . ^_ 
iDavis et al. . 2008), gene flow (ISultan and Spencerl . 120021) . and local adaptation 
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(e.g. lHanski and Gilpinl . Il998t Uoshi et all 120011) . Much interest has focused 



on the effect of the spatial structure of mctapopulations and how local popu- 
lations are connected through dispersal. Connectivity among subpopulations is 
also increasingly emphasized in management and conserv ation, e.g. to prevent 
fragmentation of landscapes ( Crooks and Saniavanl 2006) and in the design of 
protected areas and na ture reserves ( van Teeffelen et al.r2006 ). 

Early models (e.g. Levins! . 19691 Hanski . 1999 ) assumed identical disper- 
sal probability among habitat patches. The initial focus of spatially explicit 
metapopulation theory was t o explore processes that generate sp atial patterns 



in homogeneous landscapes (jHanskil . l2002t iMalchow et all 120081 ). Later, spa- 



tially explicit models were designed to let dispersal proba bility be a function o f 
patch size or the distance between local habitat patches (IHanskil . 1 19941 . 12002 ). 



One aspect of dispersal that only has been implicitly included in realistic mod- 
els but not studied in isolation is when dispersal is asymmetric. Asymmetric 
dispersal is expected for many metapopulati ons, e.g. where disp ersal is domi- 



nated by wind transport of pollen and seeds ( Nathan et all 1 2 lh T and for ma- 



rine s pecies with spores and larvae transported by ocean currents (jWares et al 
20011) . Consequently, it is important to understand how asymmetric dispersal 



may affect the dynamics and persistence of metapopulations with potential im- 
plications for the desig n of nature reserves. Some stu d ies have considered asym- 
metric dispersal (e.g. Pulliam and DanielsonL 1 1 9 9 lL iKawecki and Holt . 2002 : 
Artzv-Randrup and Stonej . 2010l) but have not analysed effects on metapop- 



ulation viability. 

In a recent contribution a conceptual model was develop ed to explore the ef- 

fects of dispersal asymmetry on metapopulation persistence (jVuilleumier and Possingham , 
20061) . The viability of metapopulations was investigated for different dispersal 



patterns randomly connecting pairs of patches through eithe r unid irectional or 
bidirectional dispersal routes. IVuilleumier and Possinghaml ([20061 ) concluded, 
that asymmetric dispersal leads to a distinct increase in the extinction risk of 
metapopulations. In a similar study iBode et al. ( 2008 ) investigated correlations 
between metapopulation viability and statistics of the dispersal network; they 
also found that asymmetric dispersal links resulted in higher extinction risk. 
Another very recently published work investigates metapopulation viability for 
a selection of asymmetric dispersal patterns and supports the findings of previ- 



ous works ( Vuilleumier et al. . 2010l) 



The main objective with this study is to isolate the effect of dispersal asym- 
metry from other properties of the metapopulation network. When chang- 
ing the degree of symmetry of dispersal networks this generally may simul- 
taneously influence the number of isolated patches and other aspects of the 
network suc h as the balance of dispersa l in th e individual patches (see e.g. 



Figure 4 in Vuilleumier and Possingham , 20061 ). Since metapopulations are 
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known to be sensitive in particular to the density of the dispersa l network 
dBarabasi and Oltvai , 2004 ) , the existence of closed cycles of d ispersal (lArmsworth , 



2008; 



20021). and the hierarchy of disp ersal in directed networks (IBode et al 
Artzv-Randrup and Stonel , 120101) these secondary implications could confound 
any effect of asymmetric dispersal. We resolve the problem by restricting our 
main analysis to regular networks. 

In this paper we in particular analyse the effect of asymmetric dispersal 
on metapopulation persistence in more detail, with an initial focus on regu- 
lar dispersal networks. We employ models of synthetic dispersal patterns and 
demonstrate that asymmetric dispersal per se may not lead to an increase in 
metapopulation extinction risk. The significance of our results, their conse- 
quence for general dispersal patterns and their relations to previous works are 
addressed in detail in Section [4] 



2. Material and Methods 

For ease of discussion we focus on the metapopulation model used in previous 



appro aches (jVuilleumier and Possinghamll2006tlBode et al. I . I2008M Vuilleumier et al 



20101) . This stochastic patch occupancy model connects a number of A patches 



through a complex dispersal matrix; the model is detailed in Section ^. II Within 
the scope of this work the viability of metapopulations exposed to dispersal pat- 
terns with different degree of symmetry is investigated. A consistent definition 
of the degree of symmetry and details on the dispersal patterns are provided in 
Sections 12.21 and [ 



2.1. Metapopulation model 

We consider a metapopulation consisting of N patches of equal quality, 
where, at a given time, each patch is either empty (0) or populated (1). In- 
teractions of the patches are specified by means of the N x N connectivity 
matrix D, where the elements dij 6 {0,1} determine whether patch j is con- 
nected to patch i (d^ = 1) or not (dy = 0). For ease of discussion we require 
da = for all i implying that patches are not connected with themselves. 

Building on previous works we used a stochastic discrete time model for a 
metapopulation of A^ patches an d tested metapopulation viability with respect 
to different connectivity matri ces ( Vuilleumier and Possingham , 2006t Bode et all 



2008t IVuilleumier et allboioh . The model, which is attractive in its simplicity, 



implements dispersal through the connectivity matrix D. Initially all A~ patches 
are populated. At each time step two events occur in succession: first, populated 
patches go extinct at per patch probability e. Subsequently, empty patches can 
be colonised with probability c by each incoming dispersal connection from a 
populated patch. Newly populated patches cannot give rise to colonisation of 
other patches at the same time step they have been colonised. 

In order to estimate the extinction risk of metapopulations the model is iter- 
ated T times. If any populated patch is left after th e T^ n iteration, the metapop 



ulation is termed viable and extinct otherwise. As Vuilleumier and Possingham 
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(2006) we chose the parameters e = 0.5 and T = 1, 000, and discuss the prob- 



ability of extinction of metapopulations consisting of N = 100 patches as a 
function of the colonisation probability c. 

2.2. Symmetry of dispersal patterns 

For characterisation of the symmetry properties of dispersal patterns the 
connectivity matrix D is divided into its symmetric and anti-symmetric contri- 
butions, S and A, by defining the matrix elements 

s i:j := mia(dij,dji) (la) 
dij := dij - Sij . (lb) 

Based on these matrices the degree of symmetry 7 of dispersal patterns is defined 
as the ratio of symmetric connections among all connections: 



Note that 1 — 7 is related to the asymmetry Z discussed in (IBode et all 120081) . 

By means of Equation @, the symmetry properties of dispersal patterns 
are put on a firm footing: Dispersal patterns are called symmetric if 7 = 1 and 
anti-symmetric if 7 = 0. Generally connectivity matrices D with intermediate 7 
are neither symmetric nor anti-symmetric. We term them asymmetric if 7 < 1 
corresponding to dispersal directed at least to some degree. 

2.3. Viability of metapopulations connected through regular dispersal patterns 

Previous works demonstrated that changes in the symmetry of dispersal 
patterns in particular affect the local symmetry of migrant flow, since asym- 
metry can result in donor- and recipient-domina t ed pa tches not present in 



symmetric networks (|Vuilleumier and Possinghaml 120061 ) . In order to isolate 



the effect of the degree of symmetry from these secondary effects, we focus on 
a specific set of dispersal patterns: we restrict our main analysis to dispersal 
patterns with the number of dispersal connections, L, being an integer multi- 
ple of N randomly distributed on the patches under the constraint, that each 
patch obtains exactly L/N in- and outgoing connections with defined degree of 
symmetry. The random patterns considered, hence , are regular with the con- 



nections evenly distributed t o all p atches available (jArtzv-Randrup and Stond . 



20ld ; Brandes and Erlebach , 2005 ). An algorithm efficiently generating regu- 



lar random dispersal patterns for small and intermediate L/N and arbitrary 
degrees of symmetry (7) is detailed in Appendix A Examples of random con- 



nectivity matrices generated for N = 8 and different combinations of L/N and 7 
are exhibited in Fig. [TJ Please regard that for the simulations metapopulations 
consisting of N = 100 are used resulting in connectivity matrices of dimension 
100 x 100 instead. 

The regular dispersal patterns we use here restrict our analysis to metapop- 
ulations with all patches connected at a fixed density independent of the choice 
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Figure 1: Examples of connectivity matrices D generated by the algorithm described in Section 
12.31 and in | Appendix A] for a reduced number of patches, N = 8, and different combinations 
of L/N and 7. Only non-zero entries are printed explicitly. For reasons of clarity symmet- 
ric connections are denoted by 'S' and asymmetric connections by 'A'. The colours indicate 
separated closed cycles of dispersal that can be identified in the matrices. While the connec- 
tivity matrices with L/N = 1 (upper row) are degenerate into 2 (7 = 0.0), 3 (7 = 0.5), and 
4 (7 = 1.0) clusters respectively, the clusters of all three matrices generated with L/N = 2 
(lower row) already extend to the entire metapopulation. fn spite of the fact that the matrices 
displayed here are only examples of randomly generated matrices, this trend is representative. 
For instance all simulations performed for N = 100 and L/N > 2 resulted in dispersal matri- 
ces with a single cluster only. Note that our results are based on much larger metapopulations 
consisting of N = 100 patches. 



of 7. For L/N > 2 the largest cluster extends to the entire metapopulation 
independent from t he degree of dispersal symmetry r esulting in irreducible con- 
nectivity matrices (ICaswelll . l200l[ iBode et all l2006h . For a detailed discussion 



of the impact of regularity on our results we refer to Section [ 

The viability of mctapopulations exposed to these dispersal patterns was 
tested in the following manner: a sample of 100 dispersal patterns connecting 
the N = 100 patches was generated for each combination of L/N = 1, . . . , 10 and 
10 different values of 7. For any of these patterns the viability of 10 independent 
realisations of metapopulations was tested for different values of the colonisation 
probability c according to the procedure outlined in Section 12.11 resulting in a 
statistics for a total of 1, 000 simulations on 100 randomly generated connectivity 
matrices for every choice of L/N, 7, and c. For our main analysis we record 
the number of viable metapopulations out of the 1, 000 simulations and prepare 
the results for graphical analysis. The sensitivity of this test procedure and its 
interpretation with respect to the statistics of extinction times is discussed in 
Section O 



3. Results 

For each scenario (L/N,j,c) a total of 1,000 simulations were performed. 
For straightforward statistical evaluation of the viability of metapopulations 
exposed to the respective conditions the simulation results were divided into 
three different groups, which are colour coded in the graphical presentation of 
the results: if all 1,000 simulated metapopulations either went extinct or were 
viable the scenario is coloured red or green, respectively. Otherwise, i.e. if the 
number of extinct simulations out of 1,000 is greater than but smaller than 
1,000, the scenario was coloured yellow. 

The results are illustrated in Fig. [2] The three panels in the upper row show 
the viability of the metapopulation as a function of the number of dispersal 
connections per patch, L/N, and the colonisation probability c for different 
values of 7: anti-symmetric dispersal (7 = 0.0), asymmetric dispersal with 
intermediate degree of symmetry (7 = 0.5), and symmetric dispersal (7 = 1.0). 
The lower panels of Fig. [2] contain the same results, but now analysed with 
respect to the effect of the degree of symmetry, 7, for three different values of 
L/N. In fact, for L/N > 3 no statistically significant impact of symmetry is 
observed. 



4. Discussion 

4-1. Interpretation and significance of results 

First of all the results depicted in Fig. [2] suggest that the impact of the 
degree of symmetry on metapopulation viability decreases with increasing L/N. 
Already at L/N > 3 no statistical significant impact of the degree of symmetry, 
i.e. no systematic differences depending on the degree of symmetry, can be 
detected on the basis of the scenarios and the statistical evaluation chosen. 
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Figure 2: Results on the viability of metapopulations exposed to dispersal patterns with regu- 
lar dispersal randomly generated by the algorithm described in Section l2.3l and |Appcndix A| 
In the upper row the viability is plotted as a function of the effective number of connections 
per patch, L/N, and the colonisation probability c. At every combination of L/N and c the 
viability of 1, 000 different dispersal patterns has been investigated. Green and red squares 
indicate parameters, where either all 1,000 patterns either were viable or not. The interme- 
diate region where some of the patterns were viable and others were not is coloured yellow. 
The three panels present the results for different degrees of symmetry, increasing from 7 = 
(anti-symmetric dispersal patterns) on the left to 7 = 1 (symmetric patterns) on the right 
hand side. In the lower row the simulation results are presented accordingly as a function of 
the symmetry 7 and the colonisation probability c for three different number of connections 
per patch, L/N = 2, 3 and 5. Only a vanishing impact of symmetry is observed for L/N > 3. 
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At a small number of dispersal connections per patch (L/N = 1, 2) metapop- 
ulation viability is significantly reduced for more symmetric dispersal (Fig. [2] 
lower panels). The reason for this effect straightforwardly can be understood 
from considerations concerning the structure of the underlying dispersal pat- 
terns: Let us first focus on patterns with L/N = 1. In this case a metapop- 
ulation with a symmetric dispersal pattern necessarily consists of a number of 
patches only pairwiscly connected through dispersal (Figure []}. The largest 
closed dispersal cycle (synony mous to the giant component of the dispersal net- 
work ( Berchenko et al.l . 120091) ). hence, involves only two patches. For the partic- 



ular metapopulation model applied a lower bound for the extinction probability 
of a pair of patches per time step is e 2 . On the contrary the mean size of the 
largest closed dispersal cycle estimated from the 100 dispersal patterns gener- 
ated for the same conditions but antisymmetric dispersal (7 = 0) was 62.7. 
For L/N = 2 the mean size of the largest cycles was 77.5 for the symmet- 
ric dispersal matrices generated, whereas for the asymmetric case all dispersal 
matrices already extended to the entire metapopulation (i.e. their mean size 
was 100). Hence we are faced with a percolation problem on random graphs 
I Callaway et al .. 2000), where the percolation threshold depends on the symme- 



try properties of the dispersal pattern. Analysis of the eigenvalues of associated 
state transition matrices reveals, that the mean time to extinction of a set of 
patches participating in a closed cycle of dispersal increases with the size of 
the cycle. For this reason differences in viability at small L/N are attributed 
to hierarchical differences of the generated matrices at only a few number of 
connections, name ly L/N < 3. This density is much sm aller then relevant cases 



discussed e.g. in (jVuilleumier and PossinghamL 120061 ) as will be discussed in 



more detail in Section! 

How meaningful is the statistical evaluation of the results with respect to 
the effect of the symmetry of dispersal patterns on expected extinction times 
of metapopulations? In order to approach this question we aim to derive lower 
and upper bounds for extinction times in the red and green regions of the 
figures, which then help to evaluate the graphical presentation of the results 
in more detail. If we disregard the initial time period of relaxation of the 
metapopulation to a quasistationary state, we can assume that the statistics 
of extinction times is exponentially distributed. This exponential distribution 
complies with a constant risk of metapopulation extinction per time step, which 
we call r. The chance, that a metapopulation has not gone extinct after T time 
steps then is (1 — r) T . For every combination of parameters we perform M 
simulations with M = 1, 000 in our cascQ. It is then straightforward to calculate 
the probability P(M\r) that all M simulations are viable, 

P(M\r) = (1 - r) MT . (3) 



1 For reasons of clarity we here assume that simulations are independent of one another 
although in each case 10 of them share the same dispersal patterns. This assumption, however, 
is not expected to be too extensive as the investigation of the replicate statistics at the end 
of Section |4. II suggests. 
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Figure 3: Extinction statistics for the metapopulations with different vaiues of the colonisa- 
tion probability c connected through dispersal matrices with 7 = 0.5 and L/N = 4. The 
individual lines indicate the number of non-extinct simulations (out of 1,000) as a function 
of the simulation time. The dashed line corresponds to the upper bound for the expecta- 
tion value of the number of extinct simulation for cases where all simulations went extinct, 
1, 000 exp(— 6.9x 10 — 3 i), as derived in the manuscript text. From the figure it becomes evident, 
that the number of non-extinct simulations after an initial relaxation phase indeed decreases 
exponentially in time (i.e. linear in this logarithmic plot). The upper bound approaches 1/M 
with t —> T, which is a general result for sufficiently large M and T as a Taylor expansion of 
expression {7} shows. For this reason the boundary line indeed exhibits the border between 
the cases marked red and yellow in Figure [2] 



Accordingly the chance that a simulation goes extinct during the T simulation 
steps is 1 — (1 — r) T , resulting in the probability P(0\r) of observing viable 
simulations of 

P(0\r) = (l-(l-rff . (4) 

More interesting, however, would be the expressions P(r\M) and P(r\0), the 
probability distributions of the mctapopulation extinction risk r given the fact 
that either all or none of the simulations are viable. These expressions straight- 
forwardly can be calculated using Bayes' theorem. Using uniform prior distri- 
butions we obtain 

-1 



P(r\M) = 



P{r\0) 



r\MT 



dr'{l - r') 



f dr'{l-{l- 
Jo 



MT 



and 



M 



(l-(l-rf) 



M 



(5) 



(6) 



Using a maximum likelihood approach confidence intervals for r can be calcu- 
lated. Applying a confidence level of 95% the upper bound for r in cases where 
all simulations are viable is 5.1 x 10~ 8 . As a lower bound for r for cases where 
all simulations went extinct we obtain 0.057. Since the latter result strongly de- 
pends on the prior distribution we instead use the inflection point of the sigmoid 
function (J6j) at 

1 - [(T- 1)/(MT- 1)] 1/T (7) 

as a more conservative estimate, which for the case of our simulations is located 
at approximately 6.9 x 10~ 3 . The inverse of r corresponds to the mean time 
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Figure 4: Extreme example of the variation in the number of viable replicates between the 
different samples (here: 7 = 0.3, L/N = 10, c = 0.15). In particular sample 98 deviates 
strongly from the general mean. Since we can assume that the main source of variations is 
the stochastic simulation procedure rather than qualitative differences between the random 
dispersal patterns relevant for the present study, we do not investigate the within-sample 
variations further within the scope of this work. 

to extinction. From our considerations we, hence, expect the mean time to 
extinction for the scenarios marked by red squares in Figure [2] to be below 
6.9 -1 x 10 3 w 145 and the respective value for the conditions marked green to 
be in the order of 2 x 10 7 or larger. Intermediate values are expected for the 
conditions marked yellow in the individual plots. Figure [3] demonstrates, that 
assumptions we needed to make seem to hold and that the estimates indeed 
reflect the underlying extinction statistics to a great extent. 

Obviously the classification of the conditions by the three scenarios to a 
meaningful extent reflects the extinction risks of the metapopulation in a sense, 
that Figure [3] succeeds to highlight the main results. From the bounds for the 
mean extinction times to extinction derived above for the respective classes we 
can conclude that metapopulations in the red regions almost surely go extinct 
within a short time, whereas metapopulations in the green regions are likely to 
be persistent. The yellow region decreases in range with increasing L/N. That 
is, the transition between threatened and persistent metapopulation sharpens 
with increasing L/N. 

The 10 replicate simulations performed for each parameter set and each 
dispersal pattern in addition allow to investigate and to discuss the variability 
within the sample of 100 dispersal patterns. In the regions marked red and 
green by definition all samples show the same behaviour. Detailed analysis of 
the yellow regions shows only very few cases of large variability of the number of 
extinct replicates between the samples. One example of rather high variability 
is depicted in Figure [4] Overall the differences between the random dispersal 
patterns generated for each scenario do not seem to be relevant for the present 
study, which is probably due to the decision of using regular dispersal patterns. 

4-2. Impact of regularity on the results 

So far we focused on regular dispersal patterns. This approach made it possi- 
ble to investigate the impact of the degree of symmetry of connectivity matrices 
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Figure 5: Results on the viability of mctapopulations exposed to general dispersal patterns 
randomly generated by modification of the algorithm described in |Appcndix A| The analysis 
and graphical presentation of the simulation results is accordant to the procedure described in 
the caption of Figure [2] Please note that L/N now specifies the mean number of connections 
per patch, while the actual number of dispersal links now can vary between patches. 



on metapopulation viability independently from other possibly confounding ef- 
fects, which is important in order to assess the role of dispersal symmetry for 
mctapopulations. Our results on regular dispersal patterns show a remarkably 
low effect of symmetry (7) on the viability of metapopulations at intermediate 
and high density of dispersal paths, L/N. At low L/N symmetric dispersal even 
results in a slightly negative effects on the viability. How do these results relate 
to the more general case where the dispersal network is not regular? 

In order to follow up this question we repeated the simulations accordingly, 
but now without the constraint of having regular dispersal networks. Techni- 
cally this was implemented by skipping steps 4c and 4d of the pattern gener- 
ation algorithm detailed in |Appendix A[ which then controls for the desired 
degree of symmetry only. The parameter L/N now should be understood in a 
statistical sense, such that L dispersal connections randomly were distributed 
between the N patches resulting in a mean density of L/N connections per 
patch. The results are depicted in Figure [SJ Interestingly, the minor effect of 
symmetry at low density of dispersal connections now shifts to a slight advan- 
tage for mctapopulations with a symmetric dispersal pattern. From L/N > 7 
no significant differences with respect to the simulation results based on regular 
dispersal patterns (Figure [5]) are observed. 

In non-regular dispersal patterns the existence of isolated patches not par- 
ticipating in dispersal has an impact on t he effective densit y of dispersal con- 



nections in the metapopulations (see also lBode et all I2008I ). Moreover, in the 
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case of asymmetric dispersal there exist patches that either only receive or only 
provide migrants, i.e. sinks or sources, and that cannot active ly take part in 
the mctapopulation dynamics (jArtzv-Randrup and Stond . l2010l) . Since both of 



these effects are most distinct at small densities of the random dispersal net- 
works, we assume that these differences basically drive the minor differences at 
low L/N between our results on regular and the general case of random disper- 
sal. Arguments for not assigning this effects to asymmetry in dispersal but to 
examine them separately are made in Section [5j 

4-3. Relation to previous works 

In general our results suggest essentially no direct negative effect of asym- 
metric dispersal on metapopulation viability at intermediate and high densi- 
ties of the dispersal network, at least as far as the stochastic patch occupancy 
mo del applied in this work is concern ed. This is in contrast to the findings 
in ( Vuilleumier and Possingham , 20061) where it was concluded that extinction 



ris k significantly increased when dispe rsal became asymmetric. The analysis 



in (jVuilleumier and Possinghaml 120061 ) is not restricted to cases with regular 



dispersal only, although the relaxation of regular dispersal is not not sufficient 
to explain the qualitative differences in the results as shown in the previous 
section. 



The description of the random patterns investigated in ([Vuilleumier and Possingham , 



20061) does not prov ide all information necessary for a n in-depth comparison 



with our results. In ( Vuilleumier and Possingham , 20061) the number of disper 



sal connections was chosen randomly for each of the 2, 000 metapopulations. 
Additional information provided on two particular patterns suggest that the 
densities are comparable or higher than the densities we investigated in our 
study. From our results we therefore do not expect a significant impact of 
dispersal asymmetry at these den sity of connections. 



The analysis of the results in (jVuilleumier and Possinghaml 120061 ) is based 



on the number of connected patches in contrast to our analysis using the global 
mean number of connections L/N. The statistics of the number of connected 
patches seems to differ significantly between the asymmetric and the symmetric 
connectivity matrices investigated, a phenomenon we were not able to repro- 
duce. In particular the example of a symmetric random pattern with more 
than 85 connections per patch but only 96 connected patches raises questions, 
since the largest cycle of closed dispersal in non-regular connectivity matri- 
ces we generated always extended to at least 99 patches for densities above 
7 connections per patch with a strong trend towards 100 patches with in- 
creasing density. For this reason we assume, that the effects described in 
( Vuilleumier and Possingham , 20061 ) originate from differences in network topol- 



ogy between the investigated connectivity matrices rather than differences in 
dis persal asymmetry- 
Bode et al. (120081) investigated the same metapopulation model as in the 



present work in a slightly different setup (N = 10, e = 0.4, and L/N = 2.6). In- 
stead of simulating individual realisations, the probability of metapopulations 
to go extinct within 100 time steps was calculated numerically for different 
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dispersal patterns. This method restricts the analysis to rather small metapop- 
ulations of 10 patches. Extinction probabilities were calculated for metapop- 
ulations connected throug h different dispersal patterns generated by the smal l 



world algorithm (see e.g. Watts and Strogat3 , 1998 ; Kininmonth et al. . 2009f) 



initiated with a particular s ymmetric dispersal pattern (Bode, pers. commu- 
2008 ) concluded from qualitative graphical analysis of 
that asymmetry reduces persistence and exhibits a 



nication). iBode et al 
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their simulation result; 
distinct threat to metapopulations. 

The discussion of our results in Section |4~T1 relates our graphical analysis to 
the extinction probability in a certain number of time step^fl, which allows for 
a comparison of the results. From additional simulation data we received from 
Bode it seems, that the negative effect in their approach is larger than what 
we would expect from our simulation for the general, non-regular case (Section 
14.21) . Additional simulations performed for metapopulations likewise subjected 
to non-regular dispersal patterns but reduced to the size of 10 patches indicated 
a general increase in the probability of extinction but no significant impact of 
mctapopulation size on the impact of symme try. We the refore assume, that the 
differences related to symmetry observed bv lBode et al. partly are owed to the 
fact, that the patterns in their study were generated from a particular symmetric 
starting configuration of the small world algorithm and that the similarity of 
patterns to this starting configuration correlates with the symmetry properties. 

Recently anot her work was devo t ed to the effect of asymmetry on metapop- 
ulation viability ( Vuilleumier et al. . 2010() . This work aims to cover different 
aspects of asymmetry simultaneously, which makes it difficult to ascribe the 
variety of effects observed to certain properties of dispersal matrices. One con- 
figuration, however, seems to be equivalent to the simulations we performed 
for general dispersal m atrices in Section 14.21 for anti-symmetric and symmetric 



dispersal, respectively ( Vuilleumier et al. . 2010l p. 229, Fig. 2, right column) 



The results the authors obtain on these patterns are in agreement with our ob- 
servations, that the degree of symmetry of dispersal matrices has no significant 
impact on metapopulation v i ability at intermediate density of dispersal connec- 



iity a 

tions (cp. I Vuilleumier et al. . 2010l p. 213, Fig. 6, difference between the plots 



in the right column). 



2 In our point of view a correlation between the extinction probability a nd di spersal asym- 
metry is not obvious from the Figure the authors refer to {Bode et all l20Qg| , p. 205, Fig. 
3). Bode, however, kindly provided additional data on an accordant simulation, which indeed 
shows a negative impact of dispersal asymmetry on the metapopulation extinction probability 
after 100 time steps. 

3 For the parameters marked green within 100 time units extinctions probabilities below 
1— exp(— 5.1 X 10~ s X 100) ~ 5x 10 -6 are expected, for the red regions an accordant calculation 
yields probabilities above almost 0.5. 
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5. Conclusions 



Wc investigated the consequences of the symmetry of dispersal patterns on 
the viability of metapopulations. Our analyses are based on simulations of a 
stochastic patch occupancy model. 

First we define the degree of dispersal symmetry, 7, which is based on the 
symmetry of the connectivity matrix (Equation [2]) • In order to be able to 
minimise possibly confounding effects we restrict our main analysis to regular 
dispersal patterns, where asymmetry does neither affect the homogeneity of dis- 
persal nor the local balances of incoming and outgoing dispersal connections. 
For these patterns we do not see any negative effect of dispersal asymmetry. For 
the more general case of non-regular dispersal patterns minor negative effects 
of asymmetric dispersal on metapopulation viability are confirmed, but only at 
rather weak densities of dispersal (cp. Section |4~2"|) . At these densities differences 
in dispersal symmetry generally arc accompanied by other hierarchical differ- 
ences of the dispersal network. This e.g. becomes evident from a neat exam ple 



of a two patch metapopulation investigated in detail in (jBode et al.l . 120081 p. 



208, Appendix A), where dispersal asymmetry by return results in a source-sink 
problem. 

From first instance it is not self-evident whether these accompanying effects 
are the origin or a consequence of asymmetric dispersal, since their characteristic 
strongly depends on how the system of study was co nstructed and chosen. Fo r 
realistic dispersal patterns the solution proposed in (|Vuilleumier et all . l2010f) . 



namely to investigate dispersal asymmetry independent from the discussion of 
sources and sinks, however docs not seem to work out, since these effects in 
general are strongly connected to one another. These correlations in the past 
made the investigation of asymmetric dispersal highly dependent on the system 
of study, which was the main difficulty in understanding the role of dispersal 
asymmetry. In order to resolve this problem wc suggest to discuss the symmetry 
of dispersal patterns at large scales e.g. based on a definition similar to Equation 
(|2|) and the statistics of sources and sinks, the homogeneity of the dispersal 
network, and other features characterising the local flow of migrants jointly 
instead of in isolation. 

It was the aim of the present work, to clarify the role of asymmetric disper- 
sal and its impact on metapopulation viability. In contrast t o previous studies 
( Vuilleumier et al. . 201fi Vuilleumier and Possingham , .2006; B ode et al. . 20081) 
we see only weak effects of asymmetric connectivity on metapopulation extinc- 
tion, which suggests that natural populations with asymmetric dispersal may 
not per se suffer from increased extinction risks. Instead effects observed in 
sim ulations, real world data, or in the evaluation of management strategies (see 



Haight and Travia I2008T ) might be reflected more significantly by other fea- 



e.g. 

tures of complex dispersal patterns. A promising path towards a discussion 
of potentially important features is taken in the investigations of the viability 
of metapopulat ions connec t ed through a variety of different dispe rsal patterns 
as provided in ( Bode et al. . 2008 : Artzv-Randrup and Stone . 201dl ). We expect 
that eventually only a theoretical analysis of the stochastic metapopulation 
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model applied can reveal the features relevant for metapopulation viability. 
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Appendix A. Algorithm for the generation of regular dispersal pat- 
terns 

Since we intended to compare cases primarily differing in their symmetry 
properties, we focused on regular dispersal patterns with fixed number of in- 
and out-going dispersal routes for every patch. For the connectivity matrices D 
this is equivalent to the constraint that the sums over every column and every 
row arc equal, that is 

J2d ik = J2d kj =L/N (A.l) 
for any i and j. Here L is the total number of activated dispersal routes. 
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Random matrices at arbitrary degree of symmetry that are complying with 
Equation (|A.ip are generated by the following algorithm, that is repeated until 
a matrix D with L non-zero elements is obtained: 

1. Set D = 0, generate a random matrix B G [0, l] NxN , where bij arc random 
numbers drawn independently from an arbitrary distribution. For instance 
uniformly distributed random variables are suitable here. Ensure that all 
elements of B are unique. 

2. Set diagonal elements ba to 10 for all i. 

3. Calculate the desired number of symmetric connections, n s = jL 

4. Repeat until smallest element of B is larger than 1 or J2i j dij = L: 

(a) Identify row i and column j of the smallest value of B 

(b) Set dij = 1 and b tJ = 10 (*) 

(c) If J2 k dik = L/N set b ik = 10 for every k (*) 

(d) If J2 k d k 3 = L/N set b kj = 10 for every k (*) 

(e) Switch i and j 

(f) if n s > (generate symmetric connection): 

• repeat the steps marked by (*) 

• reduce n s by 2 

else: (generate asymmetric connection) 

• set b^ = 10 

5. Reject result if J2i j dij < L. 

Note that the value 10, of course, is arbitrary. Any number greater than 1 is 
suitable to ensure that the corresponding elements of D arc not selected by the 
algorithm. This algorithm randomly orders the elements of D and activates 
them step by step. It generates random connectivity matrices with given degree 
of symmetry and it is sufficiently efficient for small and intermediate L. 

The implementation of the algorithm in FORTRAN90 is straightforward 
(compilation tested with the GNU Fortran compiler gfortran v4-3.3): 

REGULAR. F90 



FRDGRAM REGULARS 'CONNECTIVITY 



PROGRAM REGULAR-CONNECTIVITY 

GENERATION OF REGULAR RANDOM DISPERSAL PATTERNS 
TESTED WITH GFORTRAN 4.3.3 

(C) 2010 BY DAVID KLEINHANS , UNIVERSITY OF GOTHENBURG . SWEDEN 
DISTRIBUTED UNDER THE CREATIVE COMMONS ATTRIBUTION 3.0 LICENSE 



IMPLICIT NONE 

INTEGER., PARAMETER: : N=100 i METAPOPULATION SIZE 

INTEGER, PARAMETER: : MAX_REJECTIONS= 1000 ! MAXIMUM NO OF R.E.IECTED MATRICES 

DOUBLE PRECISION: : GAMMA /DEGREE OF SYMMETRY 

INTEGER :: LBYN ! NO . OF CONNECTIONS PER PATCH 

INTEGER: :D(N,N) 'CONNECTIVITY MATRIX 

DOUBLE PRECISION: : RAND (N , N) 'RANDOM MATRIX USED FOR. ORDERING OF LINKS 

DOUBLE PRECISION : : REMAININGJ3YM ! NO OF REMAINING SYMMETRIC CONNECTIONS 

INTEGER :: REJECTIONS ! COUNT NUMBER OF REJECTED DISPERSAL PATTERNS 

INTEGER :: LOC ( 2 ) ! LOCATION OF THE SMALLEST ELEMENT OF RAND 

INTEGER:: I, J ! AUXILIARY VARIABLES, USED FOR LOOPS ONLY 

LOGICAL :: GRIDOK ! CHECK IF GRID COMPLIES WITH CONSTRAINTS 

! === INITIALIZE RANDOM NUMBER GENERATOR === 
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CALL RANDOIVLSEED 
! REQUEST PARAMETERS === 

WRITE ( * , " ( A , I 4 ) " ) " REGULAR DISPERSAL MATRIX FOR METAPOPULATION OP SIZE N=",N 
WRTTE( * , " ( A) " ) " PLEASE ENTER PARAMETERS : " 

WRITE( * , " (A) " ) "DEGREE OF SYMMETRY, GAMMA (DOUBLE PRECISION, >=0 AND <=1)?" 

READ( * , * ) GAMMA 

WRITE( * (P8.5)" ) GAMMA 

WRITE (*,"(A,I4,A)")" NO OF CONNECTIONS PER PATCH, LBYN (INTEGER, >0 AND <" ,& 

&(N-l)/2," )7" 
READ( * , * )LBYN 
WRTTE( * ," ( 14)" ) 

/ === GENERATE DISPERSAL PATTERN == 
GRIDOK=. FALSE . 
RE.IECTIONS=0 
DO WHILE ( . NOT . GRIDOK ) 

/ == STARTING CONFIGURATION: == 
.' ALL LINKS INACTIVE 

D=0 

/ CALCULATE NUMBER OF SYMMETRIC LINKS TO BE GENERATED 

REMAINING_SYM=NINT (GAMMA* LBYN *N ) 

/ GENERATE RANDOM NUMBER MATRIX FOR ORDERING OF LINKS 

! (EXCLUDE DIAGONAL ELEMENTS BY ASSIGNING VALUE OF 10) 

DO 1=1, N 

DO J = l ,N 

IF ( I .NE. J)THEN 

CALL RANDOMJNUMBER ( RAND ( I , J ) ) 
ELSE 

RAND ( I , J ) = 10. DO 
ENDIF 
ENDDO 
ENDDO 

/ == ADD CONNECTIONS UNTIL NO LINKS ARE AVAILABLE ANY MORE == 
DO WHILE( ( MINVAL ( RAND ) . LT . 1 ) .AND. ( COUNT (D.EQ . 1 ) . LT . LBYN * N ) ) 

•LOCATE THE SMALLEST ELEMENT OF RAND 
LOC=MINLOC (RAND ) 

'.SET RANDOM NUMBER OF THE ELEMENT TO 10 AND ACTIVATE CORRESPONDING LINK 
RAND ( LOC ( 1 ) ,LOC( 2 ) ) = 1 .DO 
D(LOC(l) ,LOC(2)) = l 

'.CHECK WHETHER NUMBER OF DESIRED INCOMING OR OUTGOING LINKS ALREADY 
'.HAS BEEN REACHED FOR THE PATCH OF FOCUS, PREVENT FURTHER LINKS IF SO 
IF (COUNT(D(LOC ( 1 ) , : ) . EQ . 1 ) . GE . LBYN ) RAND ( LOC ( 1 ) . : ) = 1 . DO 
IF (COUNT (D ( : , LOC ( 2 ) ) . EQ . 1 ) . GE . LBYN ) RAND ( : , LOC ( 2 ) ) = 1 . DO 

'.IF SYMMETRIC CONNECTIONS ARE REMAINING : MAKE THE CURRENT A SYMMETRIC ONE, 
.'ELSE ENSURE THAT THE REVERSE DIRECTION IS NOT ACTIVATED 
IF ( REMAINING.SYM . GT . )THEN 

D(LOC(2) ,LOC(l)) = l 

RAND(LOC(2) , LOC ( 1 ) ) = 1 . DO 

IF (COUNT(D(LOC ( 2 ) , : ) . EQ . 1 ) . GE . LBYN ) RAND (LOC ( 2 ) , : ) = 1 . DO 
IF (COUNT(D ( : , LOC ( 1 ) ) . EQ . 1 ) . GE. LBYN) RAND ( : , LOC ( 1 ) ) = 1 . DO 
REMAINING_SYM=REMAINING_SYM— 2 
ELSE 

RAND(LOC(2) ,LOC(1)) = 10. 
ENDIF 
ENDDO 



/CHECK WHETHER THE DESIRED NO OF LINKS HAS BEEN GENERATED 
/REJECT AND RESTART IF NOT . ACCEPT THE PATTERN OTHERWISE 
IF (COUNT(D .EQ . 1 } . EQ . LBYN*N)THEN 

GRIDOK= . TRUE . 
ELSE 

RE.JECTIONS=REJECTIONS+l 

IF (REJECTIONS . LT . MAXJIBJECTIONS )THEN 



WRITE ( 1 

ELSE 

WRITE ( a 
WRITE ( * 
STOP 
ENDIF 
ENDIF 
ENDDO 



■ (A, 14 , A) " ) "PATTERN " , REJECTIONS , &c 
REJECTED, RESTARTING GRID GENERATION 



' (A)' 



' GRID GENERATION NOT SUCCESSFULL . " 

'PLEASE TRY LOWER LBYN OR INCREASE MAX_REJECTIONS . 



/ === WRITE D TO STANDARD OUTPUT 
WRITE (* ,*) "CONNECTIVITY MATRIX D: 
DO I = 1 ,N 

WRTTE( *," (99911. 1 ) " )D( I , : ) 
ENDDO 

WRITE( * , " (A) " ) "DONE! " 



END PROGRAM REGULAR_( CONNECTIVITY 
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